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ASYMPTOTIC  TECHNIQUES  IN  IGNITION  THEORY 


Final  Report  for  Grant  DAERO  79-G-0007 


by  Manuel  Kindelan 
Amable  Linan 

E.T.S.I.Aeronauticos 
Universidad  Politecnica 
de  Madrid. 


A3 S TRACT 

An  analysis  is  given  in  Part  I  for  the  ignition  of  a 
reactive  material  by  an  electrical  wire  heated  at  constant  power. 
The  wire  temperature  history  is  described  for  small  and  large 
values  of  the  electrical  heating  rate  and  for  solid  an  c  gaseous 
reactants;  in  the  case  of  gaseous  reactants  the  effects  of  con¬ 
vection  have  also  been  taken  into  account.  The  results  of  the 
analysis  include  closed  form  expressions  for  the  ignition  time 
in  terms  of  the  physicochemical  parameters. 

An  analysis  is  given  in  Part  II  of  the  process  of  ig¬ 
nition  of  reacting  particles  in  a  shock  tube.  An  analysis  for 
large  activation  energy  of  the  ignition  reaction  shows  that  for 
values  of  the  particle  diameter  smaller  than  a  critical  value, 
ignition  does  not  take  place;  the  ignition  time  is  calculated 
for  diameters  larger  than  the  critical. 

Analysis  of  ignition  by  hot  catalytic  surfaces  and  by 
inert  hot  spots,  and  of  Lewis  number  effects  on  the  structure 
and  extinction  of  diffusion  flames  due  to  strain,  carried  out 
under  this  Grant  are  being  published  elsewhere. 
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INTRODUCTION 

This  report  describes  the  results  of  applying  the 
technique  of  large  activation  energy  asymptotics,  to  several 
ignition  problems  of  practical  interest.  This  technique, 
which  has  been  extremely  useful  in  studying  many  problems  in 
combustion,  makes  possible  the  derivation  of  closed  form  an¬ 
alytical  solutions  for  the  ignition  time  or  other  properties 
of  interest,  although  the  original  equations  describing  the 
problem  are  complicated  non-linear,  partial  differential 
equations . 

The  first  part  of  this  report  studies  the  ignition 
process  of  a  reactive  material  when  an  imbedded  electric  wire 
is  heated  by  an  electric  current  of  constant  power.  The  wire 
heats  up  the  surrounding  reactive  material  by  conduction, 
until  the  exothermic  reaction  leads  to  a  thermal  runaway  at 
a  finite  ignition  time.  Since  the  heat  conductivity  of  the 
wire  is  often  very  large  compared  to  that  of  the  reactive 
medium,  the  wire  temperature  is  assumed  to  be  uniform.  The 
character  of  the  process  is  then  determined  by  the  ratio  of 
the  specific  heats  per  unit  volume  of  the  reactive  material 
and  of  the  heating  wires;  when  this  ratio  is  of  order  unity, 
the  radius  of  the  heated  region  in  the  reactive  material 
during  the  ignition  transient  is  of  the  order  of  the  wire 
radius,  while  if  the  ratio  is  small,  the  radius  of  the  heated 
region  becomes  large  compared  with  the  wire  radius.  In  the 


first  case,  the  chemical  heat  release  is  confined  to  a  thin 


layer  adjacent  to  the  wire  surface,  and  it  can  be  calculated 
in  terms  of  the  wire  temperature  and  the  heat  loss  to  the 
surrounding  material  where  the  chemical  reaction  is  frozen; 
the  surface  temperature  and  the  ignition  time  are  given  by 
the  solution  of  an  ordinary  differential  equation,  involving 
a  single  parameter.  In  the  second  case,  which  is  typical  when  ig¬ 
niting  reacting  gases,  the  characteristic  thickness  of  the 
reaction  region  is  of  the  order  of  the  particle  radius,  and 
the  radius  of  the  heated  region  is  large;  therefore  forced 
or  free  convection  effects  can  be  important  and  have  to  be 
retained  when  analyzing  the  outer  frozen  zone  of  - he  reactive 
medium . 

The  second  part  of  this  report  studies  the  process 
of  igniting  reacting  particles  in  a  shock  tube.  The  parti¬ 
cles,  which  are  located  in  a  holder  in  the  shock  tube,  are 
dragged,  accelerated,  and  heated  by  a  shock  wave.  The  wave 
is  reflected  at  the  end  of  the  tube,  so  that  when  it  encoun¬ 
ters  again  the  particles,  these  are  decelerated  and  further 
heated  by  the  gas  behind_  the  reflected  shock.  In  this  report 
we  analyze  motion  and  heat  transport  to  the  particles  as  a 
function  of  the  parameters  of  the  shock  wave,  and  we  apply 
the  method  of  large  activation  energy  asymptotics  to  study 
the  ignition  process.  There  is  a  critical  value  of  the  diam¬ 
eter  of  the  oarticles  such  that  for  DO  ..  ignition  does  not 
‘  cnt 


—  i  al¬ 
oe  c  u  r  due  to  rapid  heat  transfer  from  the  Darticle.  For  D>D 

crit 

there  is  a  runaway  of  surface  temperature  at  a  finite  ignition 
time.  Figures  are  provided  to  compute  the  ignition  time  as  a 
function  of  the  parameters  of  the  problem. 

An  analysis  has  been  carried  out  under  this  Grant  of 
the  ignition  by  a  hot  catalytic  surface,  taking  into  account  the 
effect  of  a  nonzero  rate  for  a  catalytic  consumption  of  fuel  at 
the  plate  surface.  The  analysis  leads  to  an  integral  equation 
that  describes  the  evolution  of  the  fuel  mass  fraction  at  the 
surface.  A  note  with  the  tittle  "Ignition  by  a  Hot  Catalytic 
Surface"  summarizing  this  work  was  written  by  A.  Linan  and 
Forman  A.  Williams  of  the  University  of  California,  San  Diego. 

It  will  appear  in  SIAM  J.  of  Applied  Mathematics. 

We  have  also  studied  during  the  term  of  this  contract, 
the  process  of  ignition  of  a  reactive  material  by  an  inert  hot 
spot.  The  result  of  this  research  was  presented  at  the  7th  Col¬ 
loquium  on  Gas  dynamic  of  Explosions  and  Reactive  Systems  in 
Gotingen,  August  20-24,  and  will  appear  in  the  Series  Progress 
in  Aeronautics  and  Astronautics. 

Finally,  an  analysis  of  the  "Lewis  number  effects  on 
the  structure  and  extinction  of  diffusion  flames  due  to  strain" 
was  presented  by  A.  Lir.an  at  the  International  Conference  on  the 
Role  of  Coherent  Structures  on  Modelling  of  Turbulence  and  Mixing, 
held  at  the  IBM  Scientific  Center,  Universidad  Autonoma  de  Ma¬ 
drid,  July  25-27,  1980.  The  Proceedings  will  appear  in  the  Lec¬ 
ture  Notes  on  Physics,  Springer. 


IGNITION  OF  A  REACTIVE  MATERIAL  BY  AN  ELECTRICAL  WIRE 


Part  I  of  the  Final  Technical  Report 
for  Grant  DAERO  79-G-0007 


by  Manuel  KindelSn 
Amable  LifiSn 

E.T.S.I.AeronSuticos 
Universidad  Politecnica 
de  Madrid. 


ABSTRACT 


An  analysis  is  given  for  the  ignition  of  a  re¬ 
active  material  by  an  electrical  wire  heated  at  constant 
power.  The  wire  temperature  history  is  described  for  small 
and  large  values  of  the  electrical  heating  rate  and  for 
solid  and  gaseous  reactants;  in  the  case  of  gaseous  re¬ 
actants  the  effects  of  convection  have  also  been  taken 
into  account.  The  results  of  the  analysis  include  closed 
form  expression  for  the  ignition  time  in  terms  of  the 
physicochemical  parameters. 


1.  INTRODUCTION 


Ignition  of  exothermic  materials  by  hot  wires  has 

.  1  2 

been  extensively  studied  both  theoretically  ’  and  exp e rimer. - 
3  4 

tally,  ’  in  order  to  improve  the  design  of  squibs  and  deto¬ 
nators  . 

Large  activation  energy  asymptotics  have  been  very 
useful  in  describing  thermal  ignition  with  different  types  of 
heating  mechanisms.  These  methods  are  used  in  the  present  re¬ 
port  to  analyze  hot  wire  ignition  for  different  values  of  r  he 
parameters  of  the  problem,  including  the  cases  of  gaseous  and 
condensed  exothermic  material. 

Previous  theoretical  analyses  of  ignition  of  solid 

4. 

materials  by  wire  heating,  include  the  studies  of  Altman  and 

4 

Grant  and  Kabik  et  al.  ,  who  used  the  concept  of  a  fixed  ig¬ 
nition  temperature  to  correlate  their  experimental  results. 

Friedman'*'  carried  out  a  theoretical  analysis  based 
on  the  more  fundamental  kinetic  parameters  of  activation  en¬ 
ergy  and  frequency  factor,  rather  than  the  more  ill-defined 
ignition  temperature.  The  ignition  criterion  used  minimizes 
the  sum  of  the  time  to  reach  temperature  T  by  the  inert  solu¬ 
tion  and  the  explosion  time  at  temperature  T.  An  explicit  ex¬ 
pression  for  the  ignition  time  is  derived  from  a  simple  ap¬ 
proximate  expression  of  the  inert  temperature.  This  solution 
is  in  good  agreement  with  the  results  of  numerical  integra¬ 


tions  . 


Thermal  runaway  leading  to  ignition  is  a  strongly 
dependent  function  of  temperature.  Thus,  as  it  will  be  found 
later,  an  analysis  such  as  Friedman's,  which  correctly  de¬ 
scribes  the  inert  temperature  rise  prior  to  ignition  and,  ap 
proximately,  the  temperature  at  which  thermal  runaway  occurs 
is  usually  adequate  to  describe  the  ignition  time  to  leading 
order . 

The  following  analysis  is  an  approximate  solution 
of  the  hot  wire  ignition  problem  which  is  "exact"  from  an 
asymptotic  point  of  view.  In  addition,  it  includes  several 
ranges  of  the  parameters  which  are  not  covered  by  the  approx 
imate  analysis  of  Friedman. 

The  problem  that  we  consider  is  that  of  a  cyiindri 
cai  wire  imbedded  in  a  reactive  material.  A  constant  rate  of 
thermal  energy  is  supplied  to  the  wire,  for  instance  by  an 
electric  current,  so  that  its  temperature  and  the  temperatur 
of  the  surrounding  material  increase  with  time. 

Initially,  the  temperatures  are  low,  so  that  the 
chemical  reaction  is  frozen.  Thus,  the  temperature  profiles 
are  simply  given  by  the  solution  of  the  inert  problem.  Howev 
er,  after  a  certain  time  the  energy  input  by  the  chemical  re 
action  becomes  important  so  that  the  temperature  of  the  reac 
tive  material  is  further  increased.  This  in  turn  accelerates 
the  reaction  rate,  thus  leading  to  ignition,  represented  by 
a  precipitous  rise  in  temperature,  at  a  finite  rime.  The  pur 
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pose  of  this  report  is  to  compute  this  igr. ition  time  as  a 
function  of  the  parameters  of  the  problem. 

In  Section  (2)  the  problem  of  ignition  of  reac ting- 
sol  ids  by  wires  is  formulated  in  non -dimensional  variables, 
and  the  significance  of  the  resulting  parameters  is  discussed 
In  Section  (3)  the  inert  solution  is  analysed  for  several  pa¬ 
rameter's  ranges.  These  solutions  are  then  used  in  Section 
(.4)  to  analyze  ignition  and  in  particular  to  compute  the  ig¬ 
nition  time  for  several  cases.  The  procedures  used  in  Section 
(4)  may  be  easily  generalized  in  order  to  analyze  cases  it. 
which  the  values  of  the  parameters  are  outside  the  range  of 
validity  of  the  cases  presented  in  Section  (4). 

In  practice,  convection  is  usually  important  in  the 
case  of  fluid  reactants,  so  that  the  ideal  solutions  presente 
in  Sections  (3)  and  (4)  are  of  little  applicability  to  the 
case  of  reacting  fluids.  In  fact,  either  forced  convection  is 
present,  or  else  the  temperatures  increments  associated  to 
the  wire  are  enough  to  generate  a  free  convective  flow  which 
can  not  be  neglected.  Thus,  in  Section  (5)  the  analysis  is 
generalized  to  include  the  effect  of  forced  convection,  and 
we  indicate  how  free  convective  effects  can  be  taken  into  ac¬ 
count.  Finally,  in  Section  (6)  the  results  obtained  in  the 
report  are  summarized  and  discussed. 

2 .  formulation 

'fie  consider  the  case  in  which  a  cylindrical  wire 
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imbedded  in  a  reactive  material,  is  heated  by  circulating  an 
electric  current  of  constant  power  through  it.  The  thermal 
conductivity  of  the  wire  is  considered  to  be  large  compared 
to  the  conductivity  of  the  reactive  material,  and  therefore 
its  temperature  is  uniform  throughout  the  process. 

A  distributed,  zero  order,  exothermic  reaction  of 
the  Arrhenius  type  takes  place  in  the  reactive  material  and 
produces  ignition  at  a  finite  ignition  time. 

Assuming  constant  values  for  the  density  p,  specific 
heat  c,  and  heat  conductivity  k,  the  energy  conservation  equa¬ 
tion  in  the  reactive  material  becomes 


A  exp{  -S  '  /  (  9  + 1 )  }  +  - — |- 

3  r  *■ 


1  39  _  3j[ 
r  3  r  ”  3  t  ’ 


(  1 ) 


where  all  symbols  are  defined  in  the  nomenclature.  The  tem¬ 
perature  increment  0  is  measured  with  respect  to  the  initial 
temperature  T  ,  the  space  coordinate  r  with  respect  to  the 

wire  radius  R,  and  the  time  variable  t,  with  respect  to  the 

2 

characteristic  heat  conduction  time  (R  /a  ). 

e 

The  boundary  and  initial  conditions  are 


0  ( t , co)  =  0(  0  ,r)  =  0 

,  d9s  30 

5  =  — "  2£  Tr 


(2) 

(3) 


where  the  second  equation  represents  a  balance  between  elec¬ 
tric  energy  input,  heat  absorbed  by  the  wire  and  heat  con¬ 
ducted  to  the  reactive  material.  The  parameter 


(4) 


6  = 


3  R2  ( Pc) e 


K  T  (pc), 
e  o  l 


is  a  ratio  between  the  characteristic  heat  conduction  tine 
and  the  characteristic  heat  up  time  of  the  wire,  and  the  pa¬ 
rameter 

( Pc) 


£  = 


(pc). 


(5) 


is  the  ratio  between  the  volumetric  heat  capacity  of  the  re¬ 
active  material  and  that  of  the  wire. 

Equations  ( 1 )  —  ( 3 )  have  to  be  solved  in  order  to  ob¬ 
tain  the  temperature  history  and  the  ignition  time  as  a  func¬ 
tion  of  the  four  parameters  A,  S',  S  and  e.  The  form  of  the 
solution  depends  on  the  relative  magnitude  of  these  parame¬ 
ters.  We  first  analyze  in  Section  3  the  inert  solution  for 
different  values  of  6  and  e,  and  in  Section  4,  we  use  this 
inert  solution  to  analyze  ignition  in  the  limit  of  high  ac¬ 
tivation  energy. 

Notice  that  e  is  of  order  unity  for  ignition  of 
solids  and  liquids  by  metal  wires,  while  £<<1  for  ignition 
of  reacting  gases.  The  value  of  6  depends  on  the  power  of  the 
ignition  source  and  the  wire  radius;  6  becomes  small  for  thin 
wires . 


3.  INERT  SOLUTION 


In  this  section  we  analyze  the  process  of  heating 
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a  material  by  supplying  a  constant  rate  of  thermal  energy  in 
put  through  a  wire  of  infinite  length  and  large  conductivity 
imbedded  in  the  material. 

This  process  is  described  by  equations  ( 1 )  - ( 2 )  wit 
A  =  0  .  The  solution'’  can  be  obtained  by  means  of  the  Laplace 
transform  technique,  for  instance,  and  it  may  be  written  as 


with 


9  = 


25 


1-exp ( -u  t )  Z 


du 


C  6  ) 


Z  =  J  (ru)  |  uYo(u)-2eY1(u)  |  -YQ(ru)  |  uJ<j(u)-2eJ1(u)  |  (7) 


A=  |uJ  (u)-2eJi  (u)  |  2+|  uY<j(u)-2eY1(u)  |  2 


(3) 


In  particular,  the  surface  temperature  is  given  by 


,c  .  85s 


1 -exp ( -u  t) 


u3A 


(9) 


which  for  small  and  large  times  results  respectively  in 


9  =  6(t  -  —  t3/2  +.  .  .  )  CIO) 

3  3/rr 

G  =T4r(l^1+t-Y  +  T^r  +  ...)  (11) 

sue  u  t 

where  Y  is  Euler  constant,  Y =0.5772.  Figure  (1)  shows  the 

value  of  (4  £0  /  5  )  for  different  values  of  £. 

s 

This  inert  solution  is  valid  up  to  the  time  in 
which  the  surface  temperature  increment  is  such  that  the  re¬ 
action  term  becomes  of  the  order  of  the  heat  transport 


. 
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terms.  When  that  temperature  increment  is  reached,  an  igni¬ 
tion  stage  follows  in  which  small  temperature  increments  of 
order  (S'  produce  changes  of  order  unity  in  the  reaction 
r ate  . 

The  inert  surface  temperature  increment  required  to 
start  the  ignition  stage,  is  generally  of  the  order  of  the 
initial  temperature,  and  therefore  in  this  section  we  analyze 
the  inert  solution  for  times  such  that  8  is  of  order  unity. 
In  the  following  paragraphs  we  study  the  inert  solution  for 
different  asymptotic  values  of  6  and  e,  and  in  each  case  we 
obtain  the  characteristic  times  necessary  to  develop  tempera¬ 
ture  increments  of  the  order  of  the  initial  temperature. 

a)  5  =  0(1 )  ,  s=0(l) 

The  surface  temperature  is  given  by  (9)  and  the 
heat  up  time  is  of  the  order  of  the  heat  diffusion  time, 
t=0(l).  Figure  1  may  be  used  to  compute  0s  as  a  function  of 
6  and  e  . 

b)  5  =  0(1),  a  <  <  1 

In  this  case  the  surface  temperature  is  given  to 
first  order  by  9s  =  t5  and  therefore  t=0(l).  The  temperature 
profile  in  the  reactive  material  is  given  by  the  solution  of 
the  heat  conduction  equation’  with  specified  surface  tempera¬ 
ture. 

c)  5=0(1),  £>>1 

In  this  case  t>;>l,  and  the  solution  to  equation 


(1)  is  obtained  by  matching  a  quasi-steady  solution  close  to  the  wire 
to  an  unsteady  solution  far  from  it.  This  type  of  solution 
exists  in  all  cases  in  which  t>:>l;  that  is,  when  the  heat  con¬ 
duction  region  is  large  compared  to  the  wire  radius. 

In  the  quasi-steady  region,  in  which  the  time  deri¬ 
vative  term  of  equation  (1)  is  negligible,  the  temperature 
profile  is  given  by 


e 


2tt 


=  -  rrtr  In  r  +  0 


(12) 


where  q  is  the  rate  of  heat  loss  from  the  wire  to  the  reac¬ 
tive  material,  q=-(30/’3r)  . 

s  . 

Far  from  the  wire,  there  is  an  unsteady  region  it. 
which  the  temperature  profile  is  given  by 


exp(-  4(7- dt’  (13) 

Matching  this  solution  for  small  values  of  r  to  the  quasi¬ 
steady  solution,  provides  a  relationship  between  the  heat 
flux  from  the  wire  to  the  reactive  material  q,  and  the  sur¬ 
face  temoerature  0  ,  namely 

s 


0  =  (In  4t-Y)  +  t— — 

s  4tt  4tt 


q(tZt.’)^l,Odt,  (14) 

Eq.  (14)  would  be  used  together  with 


d0 

-r^=  5  -  £3.  (14'  ) 

dt  7T 


and  0  (0)=0,  to  calculate  the  wire  temperature  history,  when- 
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q  is  small,  so  that  ignition  occurs  at  large  times. 

The  last  term  in  Eq.  (14)  is  negligible  if  the  re¬ 
lative  changes  in  q  are  small  compared  with  unity;  then  0  (t) 
is  given  by 


d8  4  e  0 

_ _ S  _  r  _ S_ 

dt  In  4 1  —  Y 


(15) 


In  the  case  considered  here,  5=0(1),  £>>1,  the 
power  supplied  to  the  wire  is  mostly  used  in  heating  the  re¬ 
active  material;  the  term  dQ^/dt  is  negligible  in  Eq.  (141), 
so  that  q  is  approximately  constant,  q=t6/£,  and 


0 

s 


( 16  ) 


The  characteristic  heat  up  time  is  of  the  order  of  (exp(e)). 


tion 


where 


d)  5>>1,  e=0(l) 

In  this  case  there  is  a  quasi-similarity  so lu- 


3/2 


0  =  t0  ,  Cn )  +  - -  9.  (n )  + . . . 


T  =  t  5  , 


/z 


n  = 


r-l 

2/7 


(17) 


(18) 


9  (n  )  =  4  i  erfc(n) 


(19) 


9  (  n  )  =  (3-4£)i^erfc(n)  +  2(nerfcr)-  6X°  - — - — -} 

2  /? 

Therefore,  the  surface  temperature  is  given  by  Eq.  (10)  and 

the  characteristic  heat  up  time  is  of  order  (5  *). 
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e )  S>>1,  s=0(5) 

Taking  the  limit  £-*-»  in  Eq.  (S),  the  surface  tern 
perature  simplifies  to 


e 

s 


2  6  f  l-exp( -u^t ) 

et2  u3| J1(u)2+Y1(u)2| 


du 


(21) 


The  surface  temperature  in  this  case  may  be  ob¬ 
tained  from  figure  1,  where  the  value  4  £  0  /  c  is  given  as  a 
function  of  £.  The  characteristic  heat  up  time  is  of  the  or¬ 
der  of  the  heat  conduction  time,  t=0(l). 


f)  6<<1,  £=0(1) 

This  case  is  analogous  to  c)  and  the  solution 
is  similarly  obtained  by  matching  a  quasi-steady  region  to 
an  unsteady  region  far  from  the  wire.  Since  5<<1  and  £=0(1) 
ail  the  energy  supplied  by  the  wire  is  used  to  heat  the  re¬ 
active  material,  and  therefore  cyconst.  Then  Eqs.  (14)  and 
(14’)  simplify  again  to  Eq.  (16).  The  characteristic  heat  up 
time  is  of  order  (exp(6  )). 


g  )  <5<<1,  £  =0  (  6  Ir.o  )<  <1 

As  in  the  previous  case,  it  i*s  necessary  to  wait 
times  long  compared  to  the  characteristic  heat  conduction 
time,  in  order  to  develop  temperature  increments  of  the  order 
of  the  initial  temperature.  Therefore,  the  surface  tempera¬ 
ture  history  is  obtained  by  matching  a  quasi-steady  region 
in  which  the  temperature  profile  is  given  by  Eq.  (12)  to  an 


-11- 


unsteady  region  described  by  Eq.  (13);  the  wire  temperature 
is  given  by  Eqs.  (14)  and  (14').  However,  in  this  case,  the 
values  of  6  and  e  are  such  that  the  energy  used  to  heat  up 
the  wire  is  of  the  same  order  than  the  energy  used  to  heat  up 
the  reactive  material,  and  therefore  all  three  terms  in  the 
surface  energy  balance  given  by  Eq.  (3),  or  Eq.  (14'),  have 
to  be  retained. 

Because  we  have  to  wait  times  of  order  1/6  for  6 

s 

to  have  increments  of  order  unity,  it  turns  out  that  the  last 
term  in  Eq.  (14)  is  small,  of  order  -l/ln5,  relative  to  the 
remaining  terms,  and  thus  q  is  given  to  leading  order  by 

q  =  4  77  9s/ln  t  ,  (22) 

and  8^  is  given  by  Eq.  (15),  which  using  the  definitions 


T 


6  , 


(  23  ) 


simplifies  to 


therefore 


d6s 

-7-=-  +  X  0  =  1  ; 

dT  s  ’ 


8  =  1 ( 1 -exp ( 6 1 )  ) 


(24) 

(  25  ) 


The  characteristic  heat  up  time  t 


is  of  order  (6  )  . 


h)  0<  < 1  ,  e  =  0 ( 6  )  << 1 

This  case  can  be  considered  a  particular  case  of 
the  previous  one  when  X<<1  so  that  in  first  approximation  ail 
the  thermal  energy  supplied  through  the  wire  is  used  to  heat 


a 


ud  the  wire,  and  0  =6t;  when  two  terms  are  retained  in  an  ex- 

s 

pansion  of  Eq.  (25)  for  small  X,  the  surface  temperature  is 
given  by 

=  •  (25’> 

_  ^ 

The  characteristic  heat  up  time  is  of  order  5 

4.  IGNITION  STAGE 

In  this  Section  we  analyze  this  ignition  stage  for 
8 ' >> 1  ,  using  asymptotic  methods  in  the  limit  of  high  activa¬ 
tion  energy  8 ' . 

Initially,  there  is  an  inert  stage  during  which  the 
reaction  term  in  Eq.  (1)  is  negligible,  and  the  inert  solution 
analyzed  in  Section  3  is  valid.  After  this  stage  there  follows 
an  ignition  stage  in  which  the  temperature  of  the  reactive 
material  is  such  that  the  Arrhenius  reaction  term  becomes  of 
the  same  order  than  the  transport  terms.  Then  small  tempera¬ 
ture  increments  of  order  1/8'  produce  changes  of  order  unity 
in  the  reaction  term;  the  energy  released  by  the  exothermic 
reaction  further  increases  the  temperature  of  the  material 
This  self -accelerating  process  leads  to  ignition,  characteri¬ 
zed  by  a  precipitous  rise  in  temperature  at  a  finite  time. 

As  discussed  in  the  previous  Section  there  are  dif¬ 
ferent  cases  depending  on  the  relative  values  of  the  parame¬ 
ters  6  and  £.  We  will  carry  out  detailed  analyses  of  six  of 
these  cases  to  show  the  methodology  used  and  to  discuss  the 


different  structures  of  the  solutions.  Similar  procedures  may 
be  used  to  calculate  the  ignition  time  in  other  cases. 


a)  6  =  0(1)  ,  e  =  0(l) 

In  this  case  the  inert  temperature  of  the  reac¬ 
tive  material  increases  by  a  quantity  of  the  order  of  the 
initial  temperature  for  times  of  the  order  of  the  character¬ 
istic  heat  conduction  time  (see  Section  3. a).  For  those  times, 
t=0(l),  the  size  of  the  heated  region  in  the  reactive  mate¬ 
rial  is  of  the  order  of  the  wire  radius  (Figure  2). 

From  Eq.  (1)  it  is  easy  to  show  that  the  character¬ 
istic  ignition  temperature  8^,  which  makes  the  reaction  term 
of  the  same  order  than  the  transport  terms  is  given  by 

exp{ 8 ' /( 9^+1 ) } =A/5  .  (26) 

Introducing  this  definition  and  using  for  convenience  as  new 
variables 

t  =  6t  ,  x=(r-l)/o”  .  (27) 


Equations  (1)  and  (3)  become 


i£  ill 

3T'9X2  + 


1 

/6"+x 


39  f  B  ’  6  8l^i 

3x  +  expU.+:  9  + 1 ' 

1  1 


(28) 


d6s  _  2£  3_9 
dt  3x 


(  29  ) 


In  the  limit  of  high  activation  energy,  the  reaction  is  fro¬ 
zen  for  9<9^,  while  it  goes  to  infinity  for  9>91  .  To  leading 
order,  the  ignition  time  is  simply  given  by  x  ,  the  time  a; 
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which  the  inert  surface  temperature  reaches  61 .  Equation  (9) 
or  Figure  (1)  may  used  to  compute  as  a  function  of  6,  £ 
and  6  ,  which  in  turn  is  related  to  A  and  S'  through  Equation 

(  26  )  . 

For  times  smaller  than  T  there  is  an  inert  stage 
which  is  followed  by  a  short  transition  stage  close  to  x  , 
which  leads  to  a  runaway  of  surface  temperature. 

During  the  transition  stage  there  is  an  inner  thin 
reaction  none  close  to  the  wire  where  the  temperature  differs 
from  5,  by  a  small  quantity  of  order  (£  ),  where 

3=3'/(ei+l)2.  (30) 

This  region  has  to  be  matched  with  an  cuter  transient-diffu¬ 
sive  tone  whose  size  is  of  order  ( 3"  ’  2 )  . 

To  analyze  the  reaction  zone  let's  define 

=  3(9-9)  .  (31) 

and  let's  expand  the  inert  temperature,  9_,  for  times  close 
to  t.  and  for  small  values  of  x 
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(28)  and  expanding  for  large  values  of  8,  ir  is  observed  that 
the  characteristic  variables  to  analyze  the  reaction  zone 
during  the  ignition  transient  are 

c-c1  =  dTS(T-T1)  ,  c  =  -  dx  6  X  (34) 


where 

a  =  -  In ( 8  d  ^ / 2 )  .  (35) 


With  these  definitions  Equations  C28)  and  (29)  result  ir. 


2  - — 4.  exD  (  a  +  -  £  )  (36) 

3q 


p 


(  2  7  ) 


where  a  single  parameter  appears, 

d_  _  dT(Tl} 

?  2e  de  d- ( ~  j  )  - i 

Figure  (3)  shows  d T / 4 £  as  a  function  of 
of  £  . 


(33) 

for  different  values 


Equation  (36)  has  to  be  solved  with  the  boundary 
condition  (Q7)  and  with  the  condition  obtained  from  matching 
to  the  outer  transient-diffusive  zone,  namely 


1  in 


a  ib 

oT 


0 


(33) 


Integrating  (36)  once,  and  using  (33)  yields 


ii> r  =  1  -  v'l  -exp (o  +  ib-h  )  •  (40) 

Evaluating  this  expression  at  the  surface  and  using  Equation 
(37)  results  i  r. 
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P  =  1  -  /l-exp(a  +  vg)  .  (41) 

Numerical  integration  of  this  equation  provides  the  surface 

temperature  history,  and  in  particular  the  ignition  time  as 

a  function  of  P,  defined  here  as  the  time  C.  when  the  ex- 

ign 

press  ion  inside  the  square-root  becomes  negative.  During  a 

_ 

short  transient,  when  a-C.  ^6  ,  the  time  derivative  term 

ign 

which  was  neglected  in  Equation  (36)  should  be  retained  to 
describe  the  thermal  runaway  at  the  ignition  time. 

Figure  (4)  shows  the  surface  temperature  history 
for  some  values  of  ?,  and  Figure  (5)  shows  the  ignition  time 
as  a  function  of  (log?).  These  results  may  be  approximated 
by  the  equation 

a  =-  In (It. 56/?) 

which  exhibits  the  correct  asymptotic  behaviour  for  small  and 
large  P . 

Once  ijj  ( a )  is  known  from  Equation  (41),  Equation 
(4Q)  may  be  integrated  to  derive  the  temperature  profile  and 
in  particular  ’He,00),  the  apparent  increase  in  wire  tempera¬ 
ture  due  to  the  chemical  reaction. 

b  )  5 > > 1  ,  £  =  0(1) 

In  this  case,  the  characteristic  inert  time  for 
a  temperature  increase  of  order  unity  is  small,  of  order 
(5  ),  compared  to  the  characteristic  heat  conduction  time 
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(see  Section  3  .  d )  .  Thus,  the  heated  region  in  the  reactive 
material  is  small  of  order  (5  '  J  )  ,  compared  to  the  wire  ra¬ 
dius  (see  Figure  6). 


The  analysis  of  this  case  p. 


.els  that  of  the 


previous  Section.  As  before,  there  is  a  transition 


whose  structure  is  described  by  a  thin  reactive-diffusive 


lose  close  to  the  wire  with  size  or  o: 


!er  (  £  ^5  ’ ' ) ,  which 


is  matched  to  outer  transient-diffusive  zone  with  size  of 
order  (3  ’ “ C  ‘  5  )  . 

I '  e  equations  to  solve  are  (2-)  and  (29),  which 
s i mol i f v  for  d  large  because  then  the  convective  term  ass: 


iated  to  curvature  in  Eq.  (23),  and  the  energy  spe 


er.t  an 


heating  the  reactive  material  in  Eq.  (23)  are  negligible,  a 
may  be  observed  in  Figure  (6). 

Introducing  the  stretched  variables  (21)  and  (24) 
in  E q s .  (23)  and  (29)  with 


=1  ,  d  =1  , 


:x“  -11'  L  i' 


is  obtained  from  Eos.  (17)-(19),  results  in  Eqs .  (35) 


=>  •:-=c 


Integrating  (26)  once  and  using  conditions  (33)  and  (42) 


yields 


V  r  ( a  ,  0  )  =  1  -  v  1  -  e  x  p  (  a  ) 


and  therefore  the  ignition  time  may  be  taken 
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t  . 
ign 


In (  200  1/tt  ) 


(45) 


where  6  is  a  function  of  6,  A  and  >5  given  by  Eq.  (26) 

c )  6<<1,  £-0(1),  0 ' =0 ( o~ 1 ) 

This  condition  which  corresponds  to  case  (f)  of 

the  previous  Section  is  used  to  introduce  the  analyses  of 

those  cases  in  wh'‘ch  t.  >>1,  and  therefore  the  structure  of 

ign 

the  solution  consists  of  a  quasi-steady  zone  close  to  the 
wire  and  an  unsteady  zone  far  fro...  it. 

Contrary  to  the  previous  cases,  the  reaction  zone 
is  not  thin  compared  to  the  wire  radius.  However,  since  the 
unsteady  effects  are  only  important  far  from  the  wire,  where 
the  temperatures  are  low  and  therefore  the  reaction  frozen, 
the  chemical  reaction  is  confined  to  a  thick  quasi-steady  re¬ 
action  zone  . 

To  analyze  this  case,  let's  introduce 

6  =  - - - j  >  4>  =  SC9-9s),  K  =  A  3  exp  (-  ("A ^~H~)  *•  U6  ^ 

(  8  + 1 ;  s 

s 

where  0  is  of  order  unity  in  the  reaction  zone.  Thus,  Eq.  (1) 
in  the  limit  of  high  activation  energy  simplifies  to 

— — +  KexD(4>)=0  (A"7) 

3r2  r  9r 

which  replaces  in  this  case  the  Eq.  (36)  describing  the  re¬ 
action  zone  structure  in  the  previous  case.  To  solve  this 
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although  the  terms  In4-Y  in  the  denominator  of  Eq.  (52)  are 
not  always  consistent  with  the  order  of  approximation  used. 
The  boundary  condition  at  the  surface  (Eq.  (3))  results  in 


94> 

(D  =  3— 
s  3r 


a  de 

(-jr-4). 


u  =  K 


(  53  ) 


From  (50)  and  (53)  an  equation  is  derived  to  describe  the 
surface  temperature  evolution,  namely 


2  ±  vT 


2+w„) “-2K 


(  54) 


or 


3' 


4£(0  +1) 
s 


=  1- 


d6 

_ S')  _ 

dt  1 


n  2 


1  - 


(0  +1)~  ln4t-YJ 
s 


(6  +1) 

s 


ex 


?  (- 


1  +  0 


-1  (55) 


This  equation  describes  the  surface  temperature  evolution  i 
all  cases  in  which  ignition  occurs  for  times  large  compared 
to  the  heat  conduction  time,  and  therefore  the  structure  of 
the  heated  region  in  the  reactive  material  is  described  by 
quasi-steady  reaction  zone  near  the  wire,  and  an  unsteady 
frozen  region  far  from  it. 

_  1 

In  the  particular  case  6<<1,  £=0(1),  $’=0(o  ), 

the  heat  absorbed  by  the  wire  is  negligible  compared  to  the 
heat  supplied  to  the  reactive  material,  while  the  chemical 
heat  release  is  'of  order  0  compared  with  the  external  heat 
input,  when  the  inert  reaches  a  value  9^  given  by 
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exo 


=  JL 

'•1  +  0  J  ?  • 


(  56) 


;me  necessary  to 


From  Section  (3.f),  the  characteristic  ignition 
reach  31  is  of  order  (exp(5  ^ 1 ,  namelv 


tl  =  4  ex?(Y 


(57.) 


Introducing  in  Sc.  (55) 


with 


T  =  t/t 


1  5 


6  _  +  0  /  S 

I  s 


(53) 


i  -r  =  5  .  +  ~r  I  n 

±  S  1  4  c. 


(  59  ) 


and  taking  me 
.--1 


limit  of  high  activation  energy  S,  of  order 
(t  i ) ,  so  that  the  characteristic  tine  in  the  transition 
stage  is  of  the  same  order  than  the  characteristic  heat  up 
time,  results  in 

b 


$=Qt  exc ( $ ) 


(60) 


where 


b  = 


4E 


Q  = 


„  o  „  2 

6.85 

i _ 

4b( b-1 ) 


(  51  ) 


This  Equation  represents  a  balance  between  the  heat  released 

by  the  reaction  and  the  increment  in  heat  flux  to  the  outer 

region  as  a  result  of  the  surface  temperature  increment  i . 

The  heat  flux  from  the  wire,  lo  ,  is  constant  and  balanced  by 

s’ 

the  heat  flux  to  the  outer  region  associated  to  the  surface 
temperature  0  „ . 


i 
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The  parameter  Q  is  always  positive,  since  b  has  to 
be  greater  than  one  in  order  to  have  the  reaction  frozen  far 
from  the  wire.  In  fact,  from  (49) 


<P  =  ai  In  r  +  cons  t . 

CO 


(52) 


and  since  u  is  assumed  to  go  to  zero  for  r-*-°°  in  order  to 
write  Eq.  (51),  it  is  necessary  that  in  the  limit  r-*-°° , 

2  2  u) 

u  =  K  r  exp(<J>)=limKr  r  “+0  (63) 

and  therefore  cu^c-2.  For  t  large  of  order  t .  ,  Equations  (52) 
and  (57)  show  that 


Li 


CO 


(  54) 


and  therefore  b>l. 

Equation  (50)  describes  the  evolution  of  the  sur¬ 
face  temperature  as  a  function  of  t,  and  the  parameters  Q 
and  b.  As  t  increases,  $  increases,  until  it  reaches  a  cri¬ 
tical  value  $=1  after  which  the  reaction  releases  heat  at  a 
rate  too  fast  to  be  transmitted  to  the  exterior,  and  there¬ 
fore  Eq.  (60)  has  no  solution.  Thus,  the  ignition  time  is 
identified  as  the  time  at  which  $=1,  so  that 

Tig„MQerl/b  ■  <65) 

Through  Eqs.  (56)  and  (61)  may  be  related  to  the  para¬ 

meters  A,  S'  ,  <5  and  £. 


d)  6  <  < 1  ,  e  =  0(l) 

This  case  is  analogous  to  the  preceding  one,  but 
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since  the  characteristic  time  of  the  transient  ignition  stage 
is  short  compared  to  the  heat  up  tine,  the  inert  solution  may 
be  e  x  d  a  r. ded  for  times  close  to  t ^  ,  the  characteristic  time  to 
reach  the  ignition  temperature  c  ^  given  by  E  q  .  (56).  Thus, 


e  =  e.  +  (T-i)  +  $/s  , 

S  1  4  c 


(66) 


where  T  is  defined  in  Ec.  (58). 

Introducing  this  relationship  in  the  Arrhenius  ex¬ 
ponent  and  expanding  for  large  values  cf  3',  it  is  seen  the’ 

_  -1 

the  characteristic  chemical  time  is  cf  crier  (it)  * .  T sing 
as  time  variable 


■3  =  (  T  -  1  ) 


(67) 


and  taking  the  limit  3 in  Eq.  (55)  results 


where 


J  =  Q  exp( v+t ) 

2 

Q  =  4  0j  £ 


Thus,  the  ignition  time  is  simply  given  by 


(53) 


(63) 


-  1  -  v-r  (  1  +  1  n  Q  ) 
ign  Bo 


(73) 


which  coincides  with  oc.  (55)  in  the  limit  b-*00  (St-4-00). 


e  )  5<<1,  e  =  0  (  6  1  n 5  )  ,  B  >  >  1  n  <5 

This  case  correspond  to  the  inert  solution  pres¬ 
ented  in  Section  )3.g)  in  which  the  parameter  a  is  of  order 
unity.  The  nondimen  si onal  activation  energy  o'  is  considered 
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to  be  large  compared  to  (lno).  Therefore,  the  reaction  zone 
is  small  compared  to  the  quasi-steady  zone,  whose  size  is  of 
order  ( ln6  )  . 

Defining  the  characteristic  ignition  temperature 

61  ,  as 

expljl^-)  =  A  Un  d)2/2Sei2  (71) 

1  1 

which  makes  the  diffusion  and  reaction  terms  of  the  same  or¬ 
der,  and  introducing  (58)  in  the  reaction  term,  with 

9  =6+6exD(-At,6)  ( t  - 1 ,  )  (72) 

I  s  1  i  l 

which  is  obtained  by  expanding  (25)  close  to  t^,  it  is  found 
that  the  characteristic  variables  in  this  case  are  $  ,  as  de¬ 
fined  in  (58)  and 

0  =  S<5  exp( -A6  t^  )  (t-t,)  .  (73) 

The  characteristic  heat  up  time  is  of  order- (6  1 )  ,  and  there 
fore  the  length  of  the  transition  stage  to  ignition  is  of  or 
der  (8  ^ ) . 

Introducing  these  variables  in  Eq.  (55)  and  taking 
the  limit  6-»-0  ,  6'=0(ln<5),  £=0(<5ln6),  results  in 

/I  -exp(  $  +  c  )'  =  1+  ||  (l-  jr-)  (74) 

whose  solution  yields  the  evolution  of  the  surface  tempera¬ 
ture  with  time,  as  a  function  of  the  parameters  of  the  prob¬ 
lem.  It  is  observed  that  this  equation  coincides  with  Eq. 


an  d 


c 


(41)  if  the  parameter  P  is  replaced  by  (1-1/15^),  and  if  the 
new  definitions  of  a  and  $  are  used.  Therefore  figure  (  5  i  mav 
be  used  to  evaluate  the  ignition  time. 


f)  3<<1,  £=0(6) 

In  this  case  the  inert  surface  temperature  is 
described  by  E  q .  (25'),  and  the  heat  flux  to  the  outer  region 

u)^,  is  negligible  compared  to  the  heat  absorbed  by  the  wire. 
This  solution  remains  valid  up  to  the  time  in  which  the  heat 
released  by  the  reaction  equals  the  heat  transmitted  to  the 
exterior,  and  therefore  the  heat  flux  from  the  wire,  it  ^  ,  be¬ 
comes  zero.  This  time  is  considered  to  be  the  ignition  time, 
and  is  obtained  by  equating  to  zero  the  expression  inside  the 
square  root  in  Eq.  (55).  Thus 

r  „  S'  ^  T  i  2  n  >2  A  8’ 


-)  *T 


(  0  t  .  + 1  )  2  n''igr.  “  ( 5 1  .  +1) 

ign  5  ign 


-exp(-Y7^ - J  (75) 

'"ign 


The 

solution 

of  this 

equation 

•u 

M 

O 

v  i  d  e  s 

t  . 

-gn 

as  a  fun 

cticn 

the 

paramet  er 

s  A  ,  8 ' 

and  0.  A  f 

o  r> 

r  ign  ’ 

the 

reaction 

none 

ceases  to  be 

quasi-s 

teady  and 

Eq  . 

(55) 

is  no 

longer 

valid 

5 .  FORCED  CONVECTION 

The  analyses  presented  in  the  preceding  sections 


nay 

b  e 

generali 

zed 

to 

include  the  effect  of  fluid 

flow  norma 

t  o 

the 

wire,  w  h 

en 

the 

Peclet  number  is  small  and 

the  parame¬ 

ter 

s  5 

and  a  ar 

e  s 

uch 

that  increments  in  surface 

ter^  e  r  a  t  u  r  e 
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i  i.  i 


of  the  order  of  the  initial  temperature  occur  for  times  long 
compared  to  the  characteristic  heat  conduction  time. 

Let  x  be  the  space  coordinate  in  the  direction  of 
flow  measured  with  respect  to  the  wire  radius,  and  let  Pe<<! 
be  the  Peclet  number  associated  to  the  convective  velocity 
U^.  The  Ossen  form  of  energy  conservation  equation  Cl) 

A  exp{-3'/(0  +  l)}  +  7^9  =  +  P  e  v  •  V  9  (76) 

can  be  used  in  first  approximation  to  describe  convective  ef 
fects.  This  equation  has  to  be  solved  with  conditions  (2)  an 
(3)  . 

In  Section  (3)  the  inert  solution  in  absence  cf 
convection  is  analyzed,  and  the  characteristic  time,  t_,  ne¬ 
cessary  to  achieve  temperature  increments  of  order  of  the 
initial  temperature  is  obtained  as  a  function  of  5  and  c.  As 
discussed  in  that  Section,  the  structure  of  the  temperature 
profile  for  times,  t  ,  large  compared  to  the  heat  conduction 
time  (t  >>1),  is  composed  of  a  quasi-steady  region  close  to 
the  wire  and  an  unsteady  region  far  from  it  where  the  temper 

atures  are  of  order  (1/lnt  ). 

•  c 

From  E  q  .  (76)  it  is  observed  that  the  effect  cf 

convection  depends  on  the  relative  magnitude  of  Fe  and  t  . 

In  fact,  the  convective  terms  in  the  unsteady  zone  are  of 

order  (?e//t  ),  while  the  transient  and  diffusive  terms  are 

c 

- 1  -i/o 

of  order  ( t ^  ).  Thus,  for  Pe  =  0 ( t  ^  ‘)  the  convective  and 
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unsteady  teras  have  to  be  retained  in  the  unsteady  region  f 
-1/2  -1/2 

r=0(t  ).  For  Pe<<t  ,  the  effect  of  convection  is  ne 

c  c 

ligible  and  the  solutions  presented  in  Sections  (3)  and  ( u ) 

-1/2 

are  valid.  For  ?e>>t  the  unsteady  term  in  Ea.  (75)  is 

negligible  and  there  is  a  quasi-steady  solution  in  which  th 
convective  terms  are  balanced  by  the  diffusive  terms.  For 
long  times,  this  solution  approaches  a  steady  state  in  whic 
the  heat  supplied  by  the  wire  is  ccr. v acted  away  by  the  flow 

_  i  /  2 

To  analyze  the  inert  solution  for  Pe>>t  '  ,  _ e t 

J  c 

introduce  as  new  dependent  variable 


$  =  6exp(-?ex/2)  , 


(  ?7) 


and  let's  assume  that  the  time  derivative  te  vn  is  neglige! 
Thus,  Eq.  (76)  simplifies  to 

32®  1  35  Pe° 


3r 


dr 


0  =  0 


(7c) 


whose  '  solution  with  conditions  (2)  and  (3)  may  be  written 

d  9 

(6  — -— )  K  (rPe/2)ext>{Pex(r-l)/2r} 

6  ,  - —SI - 2 - - -  .  (  7  s  ) 

6  '  K1(Pe/2)  -  |  K  c  (  ?  e / 2 ) 

Thus,  the  surface  temperature  history  is  given  by 


wh  ere 


9g  =  j  ( 1 -  exp ( -At)}  , 


K  (Pe/2) 

A  =  ?®  £  (Pe/2)  " 

o 


(  SC) 


(  SI  ) 
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In  the  limit  Pe-*0, 


A  =  4  £  { ln( 1 5 / ?  e  2 )  -2Y}-1 


(  32) 


so  that  the  characteristic  time  to  reach  the  steady  solution 
is  of  order  {-ln(Pe)/E}.  If  5/ A  is  small  compared  to  unity, 
the  temperature  increment  of  the  steady  state  is  small,  and 
therefore  no  ignition  occurs. 


-  _  i  /  2 

For  ?e=Q(t^  -  )  ,  there  is  a  quasi-steady  region 


close 

to 

the  wire  in  which 

both  the  convective 

and  unsteady 

terms 

are 

negligible  ,  and 

an  unsteady  zone  far 

from  it  in 

wh  i  ch 

box 

h  effects  have  to 

be  retained. 

In  the  unsteady  zone,  Eq.  (75)  has  to  be  solved 
ana  matched  to  the  quasi-steady  solution.  Using,  for  instance 
Laplace  transforms  the  temperature  profile  in  the  unsteady 
region  is  found  to  be 


9 


0  = 


exp( Pex/2  )  f  q(t)exp{-?e  ( t-T  ) /4  }  exp {■ 


4  IT 


/  4  (  t  -  T  ) 


t-T 


Evaluating  this  expression  for  small  values  of  r,  provides  a 
matching  condition  to  the  quasi-steady  region,  namely 


rt 


c ( T ) c  xp { -Pe  (t-r)/^}-c(t)  J . 


'O 


t-T 


(  3u  ) 


For  q(t)=const. ,  wren  for  instance  d<<l,  e=0(l),  so  that  all 
the  heat  supplied  by  the  wire  is  used  to  heat  up  the  reactive 
material,  this  condition  simplifies  to 

9  =  {ln( 15/Pe2)  -  E. (?e2t/4)  -  2Y} 

S  4  T7  1 


(  35  ) 
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Eq.  (34)  or 

(35)  substitutes  Sq. 

(  14 

)  as  the  mate 

hing 

rei  a 

tionship  that  together  with  Sq. 

(  14  ' 

)  is  necessar 

y  :c 

ob  t  a 

the  surface 

temperature  history. 

For 

instance,  in 

the 

case 

:<<1  ,  e=C( 1  )  , 


=  -Ar  { 1  n  (  1 5  /  ?  e  “  ) 


-  E  ,  (  ?e  -./*)  -  2Y  } 


which  in  the  limit  (?e“t)-*0,  coincides  with  Sc.  (15),  and  in 
the  limit  (Pe^t)-*-00  coincides  with  Sc.  (80)  with  the  parame¬ 
ter  A  given  by  S q .  (32). 

The  ignition  analyses  preset,  tac  in  Sections  (  3  .  c ) 

- ( 3.5)  may  also  he  generalized  to  include  the  effect  of  con¬ 
vection.  In  fact,  Sc.  (55)  is  still  valid  to  describe  the 
surface  temperature  history,  if  the  value  a:  given  by  S q . 
(52)  is  substituted  by  the  new  matching  relationship 

2  a  9 

w«,  - - 2 - ^ - 2 -  '  (  6  7  } 

ln(16/Pe  )-Z„  (?e  t/4)-2'f 


5  .  CONCLUSIONS 

The  process  of  ignition  cf  a  reactive  material  by 
an  imbedded  electric  wire  through  which  an  electrical  currerr 
of  constant  power  is  applied,  has  been  analyzed  by  means  of 
high  activation  energy  asymptotics. 

It  is  found  that  the  temperature  history  and  there¬ 
fore  the  ignition  time  is  a  function  of  four  non-dimensional 
parameters:  a  Damk holer  Number  A,  an  activation  energy  S', 
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an  energy  input  5  and  the  ratio  of  volumetric  heat  capaci¬ 
ties  £  . 

The  problem  has  been  analyzed  in  the  limit  of  large 
3',  and  it  is  found  that  the  solution  exhibits  different  be¬ 
haviour  depending  on  the  relative  magnitudes  of  S',  0  and  Z  . 

Some  of  these  limiting  solutions  are  derived  in  Section  4, 
leading  to  closed  form  solutions  for  the  ignition  time.  The 


values 

of  5  and  £ 

cons 

idered  cover 

both  the  cases  of 

vo  0  a  f  t 1  v  e 

solids 

and  gases. 

The 

same  methods 

presented  in  that 

Section 

may  be 

used  to  cb 

tain 

other  iimiti 

n g  solutions. 

Finally 

the 

effect  cf  cc 

nvectior.  which  she 

UiC  z e 

import 

ant  in  tne 

case 

of  reacting 

gases,  cs  a r. a  1  y z e c 

in  Sec- 

t ion  ( 5  )  . 
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LIST 


A  , 
B  , 


c  , 


K  , 
k  , 


L 


?e  , 


OF  SYMBOLS 

9 

Dimensionless  frequency  factor,  zR~Qpe/TQK  . 

Thermal  energy  input  to  the  wire  per  unit  volume. 
Dimensionless  parameter  defined  in  Eq.  (61). 

Heat  capacity. 

Dimensionless  parameters  defined  in  Eq.  (33). 
Dimensionless  frequency  factor  defined  in  Eq.  (46). 
Coefficient  of  thermal  conductivity. 

Dimensionless  parameter  defined  in  Eq.  (38). 

Peclet  number . 


Dimensionless  parameter  defined  in  Eqs.  (61)  or  (53 
Heat  transfer  per  unit  area  by  conduction. 


R  ,  Radius  of  the  wire. 


t  , 

u  , 

2  , 

s  , 

S'  , 
5  , 

e  , 


Dimensionless 
Temperature  . 
Dimensionless 
Dimensionless 
Dimensionless 


radial  coordinate  r/R. 

2 

time ,  ta  /R  . 

parameter  defined  in  Eq.  (15). 
radial  coordinate  Eq.  (27). 


Frequency  factor. 

Thermal  diffusivity. 

Dimensionless  activation  energy  defined  in  Eq.  (30) 
Dimensionless  activation  energy,  E/RTq. 
Dimensionless  parameter  defined  in  Eq.  (4). 
Dimensionless  parameter  defined  in  Eq.  (5). 


n 


Similarity  variable  defined  in  Eq.  (18). 
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3  , 
X  , 

a  , 
I  , 

4  , 

?  , 


4>  , 


Dimensionless  temperature,  (T-T  )/T  . 

o  o 

Dimensionless  parameter  defined  in  Eq.  (23). 
Stretched  radial  coordinate  defined  in  Eq.  (34). 


Stretched  time 

defined  in 

Stretched  time 

defined  in 

Dimensionless 

temperature 

(60)  or  (68). 

Dimensionless 

temperature 

(45)  or  (77). 

Dimensionless 

temperature 

(31). 

Eqs.  (34),  (67)  or  (73). 

Eq.  (IS). 

increment  defined  in  Eqs 

increment  defined  in  Eqs 

increment  defined  in  Eq. 


Transformed  variable  defined  in  Eq.  (43). 


Subscripts 


e 

T 


S 

0 

1 

CO 


Characteristic  variable. 

Exterior  conditions. 

Inert  . 

Interior  conditions. 

S  ur face. 

Initial . 

Characteristic  conditions  at  ignition. 
Conditions  far  from  the  wire. 
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UNSTEADY  HEAT  _J 
CONDUCTION  REGION 
lu~a/6'5 


REACTION  LAYER 
lr~  a/V06 


ASYMPTOTIC  ANALYSIS  OF  SHOCK  TUBE  IGNITION  OF  REACTIVE  PAP.TIC 

1  .  INTRODUCTION' 

An  interesting  method  to  analyze  ignition  of  reac¬ 
tive  particles  at  high  temperature,  is  the  use  of  a  reflected 
shock  technique.  In  this  technique,  the  reactive  particles, 
which  are  placed  in  a  holder  inside  the  shock  tube,  are 
"dragged",  accelerated  and  heated  by  the  incident  shock  wave. 
This  wave  is  reflected  at  the  end  wall  of  the  tube  and  when 
it  encounters  the  particles,  these  are  decelerated  and  fur¬ 
ther  heated  through  a  convective-conductive  mechanism.  These 
high  temperatures  accelerate  exothermic  reactions  so  that 
when  the  rate  of  chemical  heat  supply  is  larger  that  the  heat 
lost  to  the  gas,  ignition  occurs  at  a  finite  ignition  time. 

-i 

Cohen  and  Decker-  used  this  technique  to  measure 
ignition  delays  of  nitrocellulose  particles  in  N ^  at  gas  tem¬ 
peratures  between  600-1200  K  and  pressures  between  .1-1.5  MPa 
They  found  that  the  temperature  dependence  of  these  delays 
changes  markedly  at  750  K. 

The  purpose  of  this  Section,  is  to  use  high  acti¬ 
vation  energy  asymptotics  to  model  ignition  delays  using  the 
reflected  shock  technique.  The  model  provides  equations  for 
the  critical  diameter  and  the  ignition  time  as  a  function  of 
the  parameters  of  the  problem. 

No  attempt  has  been  made  at  this  stage,  to  compare 


the  predictions  of  the  theory,  with  the  experimental  results 
of  Cohen  and  Decker.  However,  once  calculated  the  gas  temper¬ 
ature  and  velocity  behind  the  incident  and  reflected  shock 
waves,  it  is  easy  to  compute  the  ignition  delay  time  and 
therefore  to  estimate  the  activation  energy  of  the  exothermic 
reaction  . 


2 .  FORMULATION 

We  assume  that  the  particle  temperature  is  uniform 
and  we  neglect  reactant  consumption.  Under  these  conditions 
the  momentum  and  energy  conservation  equations,  may  be  writes 
respectively  (see  Cohen  and  Decker  ) 


dU 

dt 


4  D 


P 


|  u-u ! (u-u) 


(i) 


6ssT4  Q  A  e xd ( -E / RT  ) 
CpD  +  C 


(2) 


where: 


(3) 

(4) 


The  notation  is  the  same  as  that  used  in  reference  (1),  in 
which  a  bar  is  used  to  distinguish  the  gas  properties  from 
the  particle  properties.  To  assume  a  uniform  particle  temper¬ 
ature  will  be  correct,  as  long  as  the  characteristic  heat 

0  _ 

convection  time,  t  -  =  0  C  D  /  6  k  N  u  ,  is  large  compared  to  the  char- 
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acteristic  heat  conduction  time,  t  =pCD  /'K,  in  the  particle. 

Eqs.  (1)  and  (2)  have  to  be  solved  with  the  appro¬ 
priate  boundary  conditions  describing  the  effect  of  the  in¬ 
cident  and  reflected  shock  wave  to  obtain  the  particle  tem¬ 
perature  history. 

The  nondimensional  ratio  E  /  R  T  is  usually  large. 
Then,  if  the  reactivity  is  large  enough,  the  solution  exhi¬ 
bits  a  runaway  in  the  temperature  of  the  particle  at  a  well 
defined  ignition  time.  However  as  the  reactivity  decreases, 
the  ignition  temperature  increases  and  approaches  the  gas- 
phase  temperature  T.  There  is  a  critical  value  of  the  Dam- 
konler  number,  such  that  for  values  smaller  than  the  critica 
value,  a  stationary  solution  is  reached  in  which  the  energy 
released  by  the  chemical  reaction  is  balanced  by  the  energy 
lost  to  the  gas-phase.  For  values  of  the  Damkchler  number 
greater  than  the  critical  one,  there  is  a  well  defined  igni¬ 
tion  time  which  is  derived  in  this  Section  as  a  function  of 
the  parameters  of  the  problem. 

3 ’  CRITICAL  CONDITIONS  FOR  THERMAL  RUNAWAY 

To  derive  the  critical  Damkchler  number,  let's  con 
sider  small  increments  in  the  particle  temperature  above  the 
gas-phase  temperature 

<0  =  ~  CT-T)  , 

RT 


(5) 


and  let's  analyze  the  conditions  for  the  existence  cf  a  s  • 
tionary  solution  tc  Ecs.  (l)-(2).  In  the  limit  of  high  no¬ 
vation  energies,  and  neglecting  radiation,  Sq.  (2)  may  be 


written , 


d  <>  ,  .  b 

_r  -  $  + a  e  v 
U  l 


whan  the  reaction  term  must  be  retained  because  the  tempera¬ 
ture  is  close  tc  T.  Here 


12K 

I  =  - -  t  , 

PCD" 


5  =  qa£D^|  ex?,_  J_}  (7) 

1  2  K  R  7  “  RT 


The  Nusselt  number  has  been  taker,  as  2,  since  critical  condi¬ 
tions  correspond  to  long  times  when  the  heat  transfer  by  con¬ 
vection  is  negligible. 

The  stationary  solution  of  Ic.  (£)  is  given  b-.* 

5  =  p  exp(  -4> )  (£) 


which  is  shewn  in  Fig.  1  as  a  function  cf  o 


is  c s  6 r  v “  ' 


that  there  is  a  critical  value  cf  the  Damkohler  number  {  :  __  = 

=  l/e)  such  that  for  o  >  e  no  stationary*  solution  exists:  ir. 

c 

this  case  there  is  runaway  in  particle  temperature  at  a  fi¬ 
nite  ignition  time.  For  5  <6  there  are  two  static::  a  r-*  solu- 

c 

tions  for  each  value  of  0,  cf  which  only  the  ievrer  one  is 
stable.  From,  the  definition  of  c,  the  critical  diameter  is 
given  by  the  2  e m e n c  v  result 


u ■ ^ 1  K  RT 

pea; 


exp ( E /RT ) 


When  the  effect  of  radiation  is  retained  the  crit¬ 
ical  regime  occurs  with  temperatures  close  to  the  chemically 
frozen  equilibrium  temperature  T^.  given  by 

Tv  -  -r  .  £OD  4 

1  ~  l  -  T  -  i  r- 

1  2K  f 


In  this  case  the 


regime  is  still  given  by  Eq.  (5) 


and  the  critical  diameter  by  Eq.  (9),  is  replaced  in  uqs 


3  ,  7T, 


(7)  and  (9)  by  K { 1  +  2e aDT  /  / K }  and  T  by  T  - 


4.  IGNITION  ANALYSIS 

"or  values  of  D  of  the  order  of  D  the  transient 

c 

Temperature  history  is  described  by  Eq.  (6).  For  values  of 

D  <  <  D  a  thermal  runaway  will  occur  when  the  particle  temcer- 
c 

ature  is  smaller  then  T  (or  T^.)  by  a  quantity  large  compared 
with  RT2/E. 

To  analyze  in  this  case  the  thermal  runaway  to  ig¬ 
nition,  we  carry  out  an  asymptotic  analysis  in  the  limit  of 
high  activation  energy  of  Eqs.  (l'-(2),  and  consider  t  as 
the  time  of  thermal  runaway. 

Let's  define,  T  =T,(t  ),  as  the  tenmerature  of  the 
particle  at  time  tc,  with  the  chemical  reaction  frozen.  In 
the  limit  of  high  activation  energy,  the  chemical  heat  re¬ 
lease  is  negligible  compared  ;o  the  convection  of  heat  from 
the  hot  gases,  as  long  as  Tj.<Tc<  Thus,  for  times  smaller  than 
t  ,  the  temperature  history  is  gi  ;en  by  the  solution  of  Eqs. 


( 1 ) — ( 2 )  without  the  chemical  source  term. 


For  times  close  to  t  ,  there  follows  a  sho 
of  transition  to  ignition  in  which  the  chemical  heat 
term  becomes  important  and  has  to  be  retained,  burin 
stage,  the  temperature  of  the  particle  differs  from 
small  quantity  of  order  R  /E,  but  these  small  chan 
fice  to  produce  changes  of  order  unity  in  the  reacti 
To  study  the  ignition  stage,  let's  introdu 


<t>  = 


R  T 


CT-TI)  , 


and  let's  expand  T,  for  times  close  to  t 


T  =  T  + 
c 


dTH 


(■-'c)  + 


Introducing  this  expansion  in  Eq.  (2),  and  defining 
time  variable 


dT, 


a  =  (t-t  ) 


c  a1 


( 


t  R  T  ' 
c  c 


results  to  leading  order  in 


||=  P  e  xp  (  j  +  J  ) 


with  initial  conditions 

a  ->  -»  cf>  =  0 

where 


( 


-7- 


Inregrating  Eq.  (13),  results  in  the  following  particle  tem¬ 
perature  history 

<p  =-  In ( 1  -P  e3)  .  CIS) 

Since  t  is  defined  as  the  time  at  which  thermal  runaway  oc¬ 
curs,  Eq.  (13)  has  to  satisfy  the  condition 

0  =  0,  <J>  -*■  “  (17) 

and  therefore  P  =  1 . 


Thus,  the  ignition  time  is  defined  by  the  equat 
^  exp( -E/RT,  )  =  -jr-  .  (13) 


con 


Thereto  re 

,  to  obtain  th 

v  to  derive 

the  inert  tern 

without  the 

reaction  term 

to  deduce 

the  ignition  t 

5.  INERT  SOLUTION 


To  solve  Eq.  (1),  let's  introduce  the  nondimensionai 


variao-s 


x  = 


Re 
Re  . 


(  19  ) 


where 


Re  .  = 


OD  I  U-U . 

1  l 


_4  OD 


C  3  U  Re. 2/3  ‘ 


With  these  definitions  Eq.  (1)  results  in 


(  20) 


dx 

dx 


=  -  x  ( 


'  4 


24 


1/3  _  2/3 

x  x  Re  . 


)  . 


(  21  ) 


-9- 


Eq.  (13),  defining  the  ignition  time,  may  be  written 

,2 


(  T  -  T  )  —  —  ^  £  _ 


exp( -E/RT) 


(  27  ) 


6  K  z(2+.46  Pr  '  Re  /  ^  x' 

l 


which  together  with  Ec  .  (23)  defines  x.  and  therefore  t  as  a 

ign  ign 

function  of  the  parameters  of  the  problem. 

To  solve  Eq.  (27),  one  may  take  advantage  of  the 


fact  that  the 

activation  ene 

in 

is  la 

nge  , 

and  use  the 

ing  iterative 

procedure.  Let 

*  s 

define 

^  _  Q  A  0  D  2 

-  .  55-, 

6  K  (  T  -  T  )  ,  ?  r  ‘ 

-  -  5 

O  0  -  .  Z)  z 

r,e  . 

- 

.3 

z'  r 

"  Re/ 

H  +  -46x  > 

l 

and  let's  take 

the  natural  1 

3g 

of  Eg. 

(  27  ) 

• 

RT  +  1.*  j 


However,  since  j  is  approximately  of  orde: 
large,  one  may  neglect  In  j  ,  and  obtain 


(  29  ) 

iity,  while  f  is 


ign 


(In  f ) 


( T-T) .  R( T-T) . 
l  l 


(30) 


Figs.  2  and  3  may  than  be  used  to  derive  the  zeroth  order  ap¬ 
proximation  to  x.  and  T.  .  using  these  values  in  Ec.  (28) 

ign  ign 

a  zeroth  order  approximation  to  j  is  calculated,  which  may 
then  be  used  in  Eq.  (29)  to  derive  the  first  order  approxi¬ 
mation  to  z.  .  This  iterative  orocedure  converges  very  r  a  - 
pidly  to  the  ignition  time. 


A  further  simplification  may  be  introduced  by  con- 


sidering  that  the  initial  Reynolds  number,  Re..,  is  large 
compared  with  unity,  as  will  usually  be  the  case  in  practice. 
Using  the  solutions  to  Ecs.  (21)  and  (24)  for  large  values 
of  Re.  in  Eq.  (30),  an  explicit  zeroth  order  solution  for 
the  ignition  time  is  derived,  namely 


3 1  . 
len 


1666  .  36  5  Re  .’ 


,  ..  f  T 

^  -n  |  _ 

HT-T) 


(inf) 


R( T-T)  . 


This  equation  is  valid  as  long  as  T.  <<T.  When  T.  a  o  - 

proac'nes  T,  the  critical  solution  discussed  in  (3)  should 
be  used. 

It  should  be  pointed  out,  that  in  the  case  in  which 
ignition  occurs  in  the  reflected  shock  region,  the  ignition 


ters  the  rerlectec  shoe*  ,  and  the  time  from  this  instant  to 
thermal  runaway. 

The  first  component  of  the  ignition  time  is  calcu¬ 
lated  from  the  parameters  of  the  incident  and  reflected  shock 
waves,  distance  from  the  holder  to  the  end  wall  and  Fig.  2 
giving  the  particle  velocity  as  a  function  of  the  non -dimen¬ 
sional  time  x.  Also,  from  Fig.  3  the  temperature  increment 
in  the  particle  (T-T),  existing  at  the  instant  of  encoun* er 


with  the  reflected  shock  wave,  is  calculated  as 


r unction 


of  the  initial  temoerature  increment  (T-T)..  The  resulting 

l  55 

temoerature  increment  is  used  as  the  initial  value  in  Ecs. 


(  1 ) -  (  2 )  to  calculate  the  time  from  the  instant  the  particle 


enters  the  reflected  shock  region  until  thermal  runaway  oc¬ 
curs.  Analogously,  the  value  Re^  in  Eq.  (27)  is  calculated 
from  the  particle  velocity  at  the  time  of  entrance  in  the  re¬ 
flected  shock  region,  which  may  be  obtained  from  Fig.  2. 
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